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ABSTRACT 

Pour  crystal  plasticity  codes,  the  viscoplastic  Material  Point  Simulator  (MPS)  developed 
at  Cornell  and  the  ViscoPlastic  Self-Consistent  code  (VPSCVb)  developed  at  LANE,  and  two 
elastic-viscoplastic  codes  developed  at  Drexel  University,  were  employed  to  calculate 
deformation  textures  and  mechanical  properties  of  model  polycrystalline  specimens  by 
simulating  isochoric,  free  upsetting.  Uniaxial  compression  of  a  model  sample  with  a  starting 
random  texture  of  5000  grains  was  carried  out  at  a  constant  true  strain  rate  of  0.001/s  to  a  true 
strain  of  1.0  with  0.02  strain  increments.  Material  properties  simulated  a  face-centered  cubic 
(PCC)  alloy.  Type  304  Stainless  Steel,  and  a  hexagonal  close-packed  (HCP)  material,  unalloyed 
Ti.  Both  non-hardening  and  linear  hardening  conditions  were  investigated.  Different  strain-rate 
sensitivities  simulated  deformation  conditions  appropriate  to  ambient  and  elevated  temperature 
conditions.  All  codes  permitted  use  of  the  Taylor  homogenization  hypothesis,  resulting  in  an 
upper  bound  for  the  mechanical  properties.  All  codes  produce  essentially  identical  results  for  the 
same  input  material,  homogenization  hypothesis  and  deformation  conditions.  Eor  comparison, 
one  alternative  homogenization  hypothesis  to  model  grain  interactions  was  examined  for  each  of 
the  MPS  and  VPSCVb  codes. 
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INTRODUCTION 

Over  the  last  two  deeades,  significant  progress  in  the  field  of  crystal  plasticity  modeling 
has  led  to  the  development  of  many  codes  with  proven  success  in  simulating  various  aspects  of 
the  mechanical  behavior  of  both  FCC  and  HCP  metals.  This  proliferation  makes  it  difficult  for 
researchers  to  decide  which  code  is  best  suited  for  a  particular  application.  This  study  aims  to 
provide  an  unbiased  study  of  the  capabilities  and  limitations  of  three  different  types  of  codes  by 
examining  the  results  obtained  from  each  using  identical  input  conditions.  The  codes  employed 
are  1)  the  ViscoPlastic  Self-Consistent  code^  (VPSCVb)  developed  and  maintained  by  C.  Tome 
and  R.  Lebensohn  at  Los  Alamos  National  Laboratory,  Los  Alamos,  NM,  hereafter  referred  to  as 
V;  2)  the  Material  Point  Simulator  (MPS)  code^  developed  by  P.  R.  Dawson,  D.  E.  Boyce  and 
associates  at  Cornell  University,  hereafter  referred  to  as  C;  and  3)  two  elastic-plastic  codes,  one 
each  for  face-centered  cubic  and  hexagonal  close-packed  metals,  developed  by  S.  Kalidindi  and 
associates  at  Drexel  University^’"^,  hereafter  referred  to  collectively  as  D. 

All  of  the  codes  are  capable  of  calculating  effective  stress-effective  strain  (SS)  curves,  the 
final  orientation  of  each  grain  in  the  ensemble  (TEX),  from  which  pole  figures  can  be  calculated, 
and  the  effective  Taylor  factor  (M)  for  the  ensemble.  All  employ  a  form  of  the  Voce  hardening 
law  for  hardening  on  the  active  slip  systems  and  the  V  and  D  codes  permit  the  inclusion  of  latent 
hardening  in  cases  where  more  than  one  type  of  slip  system  is  active.  Although  all  codes  are 
capable  of  including  both  slip  and  twinning  in  the  calculations,  only  slip  was  considered  in  this 
study  to  reduce  the  number  of  permutations  of  material  variables.  A  homogenization  assumption 
specifies  the  relationship  between  the  state  of  deformation  for  individual  grains  and  the  global 
state  of  deformation  for  the  ensemble.  The  Taylor  hypothesis,  in  which  the  strain  or  strain  rate  in 
each  grain  is  the  same  as  that  for  the  ensemble,  is  common  to  all  three  codes  and  gives  an  upper 
bound  to  the  resulting  SS  curve.  In  addition  the  C  code  provides  an  option  for  the  Sachs 
hypothesis,  in  which  the  state  of  stress  in  each  grain  is  the  same  as  for  the  ensemble,  giving  a 
lower  bound  to  the  SS  curve.  Although  this  code  also  permits  selected  weighted  averages  of  the 
Taylor  and  Sachs  hypotheses,  these  were  not  investigated  in  this  study.  The  Affine  linearization 
hypothesis  of  the  V  code^  was  employed  in  this  work  for  comparison. 

MATERIAE  PARAMETERS  AND  COMPUTATIONAE  DETAILS 

Material  properties  employed  in  the  calculations  correspond  to  austenitic  stainless  steel,  a 
face-centered  cubic  (ECC)  avatar,  and  unalloyed  Ti,  a  hexagonal  close-packed  (HCP)  avatar. 
Crystal  deformation  was  confined  to  {111}<110>  slip  for  the  ECC  material  and  to  the  basal  <a>, 
first  order  prismatic  <a>  and  pyramidal  <c  +  a>  systems  in  the  HCP  material.  Eor  both  ambient 
and  high  temperature  deformation  conditions,  relative  values  of  the  critical  resolved  shear  stress 
(CRSS)  on  each  of  the  slip  systems  in  the  HCP  material  were  set  to  1.0:0. 7:3.0,  respectively. 
Deformation  by  twinning  was  not  included  in  any  of  the  calculations.  Reduced  values  of  the 
CRSSs  on  each  slip  system  were  employed,  so  that  the  CRSS  on  the  reference  slip  system  (basal 
for  hep)  was  1.0. 

Different  temperatures  were  simulated  only  by  employing  different  values  of  the  stress 
exponent  (strain-rate  sensitivity)  while  the  hardening  ratios  and  CRSSs  remained  constant.  Eor 
both  the  ECC  and  HCP  materials,  high  temperature  deformation  conditions  correspond  to  a  stress 
exponent  n  =  4  and  low  temperature,  to  a  stress  exponent  of  n  =  20  in  an  equation  of  the  form 
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where  is  the  ij*  eomponent  of  the  deviatorie  strain  rate  due  to  slip  on  the  k*  slip  system, 
m|^  is  the  Sehmid  tensor  on  that  system,  Ors  is  the  loeal  deviatorie  stress  tensor,  n  is  the  stress 


exponent,  Yo  is  a  referenee  strain  rate,  and  and  r^is  a  referenee  stress  for  slip  on  the  k™  slip 
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system,  respeetively.  The  strain  rate  due  to  slip  on  all  systems  is  obtained  by  summing  equation 
(1)  over  all  aetive  slip  systems.  While  the  V  eode  employs  this  form  of  the  rate-dependent  law, 
the  C  and  D  eodes  use  the  form  in  whieh  the  resolved  shear  stress  is  given  in  terms  of  the  strain 
rate,  using  a  strain-rate  sensitivity  exponent,  m  =  1/n. 

Both  hardening  and  non-hardening  eonditions  were  investigated  using  an  empirieal  Voee- 
type  hardening  law.  On  eaeh  aetive  slip  system  these  relationships  have  the  form 

^(Y)  =  ^o  +  (^s-^o)(l-e^“')  (2) 

where  t(y)  is  the  resolved  shear  stress  at  a  resolved  shear  strain  of  y,  Xo  and  is  are  the  eritieal 
resolved  shear  stress  (CRSS)  and  saturation  stress,  respeetively  on  the  slip  system  and  a  is  a 
hardening  parameter.  The  V  eode  employs  a  version  of  equation  (2)  that  uses  the  aeeumulated 
shear  strain  in  a  grain  instead  of  y  for  the  strain  and  permits  hardening  to  approaeh  a  non-zero 
asymptotie  value,  introdueing  a  fourth  parameter  into  the  expression'.  The  present  study  treats 
only  linear  hardening,  simulated  over  the  range  of  strain  for  whieh  ealeulations  were  performed 
by  setting  is  =  1000  to  and  the  initial  hardening  rate  0o  =  a(Ts-Xo)  =  0.2xo.  For  the  HCP  material 
the  initial  hardening  rates  on  the  aetive  slip  systems  were  seleeted  in  the  same  ratios  as  the 
CRSSs. 


The  effeetive  Taylor  faetor,  M,  for  the  ensemble  of  grains  was  ealeulated  at  eaeh 
inerement  of  effeetive  strain  for  the  V  and  D  eodes.  In  the  former  ease  M  is  given  by 


M  = 


G  8 
y  y 


(3) 


where  summation  over  repeated  indiees  is  implied  and  the  deviatorie  stress  and  strain  rate 
eomponents  refer  to  global  values.  The  SS  eurves  were  normalized  by  dividing  stresses  by  the 
CRSS  on  the  referenee  slip  system  (basal  slip  for  HCP).  Equation  (3)  and  the  use  of  1.0  as  the 
CRSS  on  the  referenee  slip  system  results  in  the  value  of  M  being  numerieally  equal  to  the  von 
Mises  effeetive  stress  at  eaeh  inerement  of  strain  for  both  FCC  and  HCP  materials.  The  D  eode 
uses  a  definition  for  M  that  replaees  individual  stress  eomponents  in  equation  (3)  with  the  CRSS 
on  the  referenee  system. 

Computations  of  the  SS  eurves  were  performed  at  effeetive  strain  intervals  of  0.02  from 
an  initial  true  strain  of  0.0  to  a  final  value  of  1.0  under  uniaxial  eompression  along  the  X3  axis  at 
a  rate  of  10'  s'  .  However,  beeause  of  the  elastie-viseoplastie  nature  of  the  D  eodes,  when  the 
instantaneous  total  strain  was  less  than  0.05  it  was  neeessary  to  perform  ealeulations  at  a  smaller 
inerement  of  strain  than  for  the  rest  of  the  eurve  to  aehieve  agreement  with  the  purely 
viseoplastie  eodes  in  the  vieinity  of  initial  yielding,  an  ensemble  of  5000  randomly  oriented 
grains  was  seleeted  as  the  starting  material. 
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RESULTS 
Face-Centered  Cubic 

The  SS  eurves  for  the  FCC  calculations  produced  by  all  eodes  are  shown  in  Figure  1.  The 
legend  identifying  data  sets  specifies  the  code  (C,  D,  V),  crystal  structure  (F  or  H),  temperature 
as  high  (H)  or  ambient  (R),  hardening  (H)  or  non-hardening  (N),  and  the  homogenization 
assumption.  The  latter  is  represented  by  U  for  the  Taylor  hypothesis,  L  for  the  Sachs  and  A  for 
Affine  linearization.  The  effeetive  stress  is  numerieally  equal  to  M  for  the  V  eode. 


Figure  1.  FCC  Upper  Bound  Effeetive  Stress-Strain  Curves  from  All  Codes 


Pole  figures  from  the  calculations  shown  in  Figure  1  are  given  in  Figures  2a  and  2b  for 
hardening  and  non-hardening,  respectively,  under  ambient  temperature  deformation  conditions. 
Orientations  of  all  grains  in  the  ensemble  at  an  effeetive  true  strain  of  1.0  were  employed  to 
construet  pole  figures  for  the  deformed  materials  using  TSL^  software  marketed  by  ED  AX.  All 
pole  figures  have  the  compression  (xs)  axis  at  the  center  of  the  figure  and  intensity  contours 
represent  multiples  of  the  random  probabilities  for  the  indicated  poles. _ 


100  110  100  110  100  110 
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Figure  2a.  FCC  Pole  Figures  for  Ambient  Temperature,  Hardening  Deformation  Condition 


Figure  2b.  FCC  Pole  Figures  for  Ambient  Temperature,  Non-  Hardening  Deformation  Condition 


Figures  3a  and  3b  show  similar  pole  figures  for  high  temperature  conditions. 
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Figure  3a.  FCC  Pole  Figures  for  High  Temperature,  Hardening  Deformation  Condition 
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Figure  3b.  FCC  Pole  Figures  for  High  Temperature,  Non-Hardening  Deformation  Condition 


Figure  4  compares  the  SS  curves  calculated  with  the  C  code  using  the  U  and  L  assumptions. 
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Figure  4.  FCC  Upper  and  Lower  Bound  Effeetive  Stress-Strain  Curves  from  C  eode 
Figure  5  shows  pole  figures  associated  with  the  lower  bound  stress-strain  curves  in  Figure  4. 


Figure  5.  FCC  Pole  Figures  for  Lower  Bound  calculations  with  C  code 


Figure  6  compares  stress-strain  curves  obtained  with  the  V  code  using  the  Taylor  (T)  model  with 
the  same  deformation  conditions  calculated  using  the  Affine  (A)  option  of  the  code. 
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Figure  6.  FCC  Effective  Stress-Strain  Curves  from  Taylor  and  Affine  options  of  V  code 
Figure  7  gives  pole  figures  obtained  using  the  Affine  option  of  the  V  code. _ 


100  110  100  110  100  110  100  110 
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High  Temperature  Ambient  Temperature 


Figure  7.  FCC  Pole  Figures  from  Affine  option  calculated  with  V  code 


Hexagonal  Close-Packed 

Figure  8  shows  SS  curves  obtained  from  all  codes  for  HCP  material. 
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Figure  8.  HCP  Upper  Bound  Effeetive  Stress-Strain  Curves  from  all  codes 
Pole  figures  for  the  stress-strain  curves  in  Figure  8  are  shown  in  Figure  9a,  9b,  10a  and  10b. 


VHRNU 


Figure  9a.  HCP  Pole  Figures  for  Ambient  Temperature,  Hardening  Deformation  Condition 
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Figure  9b.  HCP  Pole  Figures  for  Ambient  Temperature,  Non-Hardening  Deformation  Condition 
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Figure  10a.  HCP  Pole  Figures  for  High  Temperature,  Hardening  Deformation  Condition 


DHHNU 


VHHNU 


CHHNU 


Figure  10b.  HCP  Pole  Figures  for  High  Temperature,  Non-Hardening  Deformation  Condition 


von  Mises  Effective  Strain 


Figure  11.  HCP  Upper  and  Lower  Bound  Effective  Stress-Strain  Curves  from  C  code 


9 


Figure  1 1  compares  lower  bound  and  upper  bound  SS  curves  for  HCP  material. 
Pole  figures  from  the  lower  bound  calculations  are  shown  in  Figure  12. 


Hardening  Non-hardening  Hardening  Non-hardening 

High  Temperature  Ambient  Temperature 


Figure  12.  HCP  Pole  Figures  for  Lower  Bound  calculations  with  C  code 


Figure  13.  HCP  Effective  Stress-Strain  Curves  from  Taylor  and  Affine  options  of  V  code 
Pole  figures  obtained  with  the  Affine  option  are  shown  in  Figure  14. 
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High  Temperature  Ambient  Temperature 


Figure  14.  HCP  Pole  Figures  from  Affine  option  calculated  with  V  code 
DISCUSSION  AND  CONCLUSIONS 

The  ratio  of  initial  flow  stresses  for  the  polycrystalline  effective  stress-strain  curves  under 
high  and  ambient  temperature  deformation  conditions  was  nearly  equal  to  the  strain  rate  raised  to 
a  power  [(1/nL)  -  (l/nn)],  where  nn  and  nL  refer  to  the  stress  exponents  for  high  (H)  and  ambient 
(L)  temperature  deformation.  Thus  the  difference  in  initial  flow  stresses  is  accentuated  by  the  use 
of  a  relatively  slow  strain  rate  in  the  calculations.  If  the  strain  rate  employed  for  the  calculations 
were  1.0  s"',  there  would  be  virtually  no  difference  in  the  initial  flow  stresses  of  the 
polycrystalline  effective  stress-strain  curves  obtained  for  the  two  deformation  conditions.  The 
differences  that  do  exist  can  be  attributed  to  a  small  amount  of  texture  hardening. 

Hardening  exhibited  by  the  effective  stress-strain  curves  calculated  using  the  Taylor 
assumption  for  the  polycrystalline  materials  is  virtually  linear,  as  is  the  input  constitutive 
equation.  However,  the  intensity  of  hardening  is  enhanced  in  the  polycrystalline  aggregate  as  a 
result  of  texture  formation.  While  the  input  hardening  rates  are  0.2(CRSS)  on  each  active  slip 
system,  the  apparent  linear  hardening  rates  of  the  FCC  polycrystal  curves  are  0.7  and  0.6  times 
the  initial  flow  stress  for  the  ambient  and  high  temperature  deformation  conditions,  respectively. 
The  corresponding  factors  for  the  HCP  polycrystal  curves  are  0.46  and  0.17.  This  enhanced 
hardening  of  the  polycrystalline  aggregate  may  also  be  attributed  to  texture  hardening  effects. 

Texture  differences  between  upper  and  lower  bound  calculations  obtained  with  the  C 
code  can  be  attributed  to  two  characteristics  of  the  lower  bound  assumption;  1)  in  the  lower 
bound  approximation  fewer  grains  accommodate  more  deformation  in  contrast  to  the  upper 
bound  where  all  grains  experience  the  same  deformation  and  2)  since  all  grains  experience  the 
same  stress,  hardening  of  a  grain  renders  it  less  capable  deforming  further  forcing  other  grains  to 
accommodate  more  of  the  deformation.  These  effects  cause  the  texture  to  be  more  strongly 
developed  in  the  hardening  cases  for  both  high  and  low  temperature  deformation.  With  no 
hardening,  all  grains  remain  active  and  reorient  as  deformation  progresses,  while  hardening 
causes  decreased  slip  activity  in  some  grains,  forcing  the  activation  of  other  grains. 

All  three  codes  produced  excellent  agreement  among  calculated  effective  stress-strain 
curves  for  all  material  and  deformation  conditions  investigated.  Pole  figures  produced  by  all 
codes  were  essentially  identical  for  the  Taylor  homogenization  hypothesis.  The  principal 
differences  in  pole  figures  for  both  FCC  and  HCP  material  occurred  as  a  result  of  different 
strain-rate  sensitivities. 
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